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I. Introduction 

T HE Global Positioning System (GPS) has reached operational status and is 
now being used internationally as a means for accurate, world-wide, all- 
weather positioning and navigation based on pseudonoise signals transmitted by 
a constellation of 24 satellites. Pseudolites (PLs) are ground-based transmitters 
that can be configured to emit GPS-like signals for enhancing the GPS by 
providing increased accuracy, integrity, and availability. 11 - 21 

Accuracy improvement can occur because of better local geometry, as measured 
by a lower vertical dilution of precision (VDOP), which is important in aircraft 
precision approach and landing applications. Accuracy and integrity enhancement 
can also be achieved by employing a PL’s integral data link to support differential 
(DGPS) modes of operation and timely transmittal of integrity warning informa¬ 
tion. Availability is increased because a PL provides an additional ranging source 
to augment the GPS constellation. 

Although the use of PLs offers many potentially significant benefits, a number 
of technical issues must also be addressed. One is the PL signal power level and 
the associated “near-far” problem that a user receiver may experience, depending 
upon the dynamic range of signal strength encountered as the distance to a PL 
changes. Other issues include deployment requirements, signal data rate, signal 
integrity monitoring, and user antenna location and sensitivity. This chapter 
addresses PLs from the perspectives of signal design considerations, integrated 
DGPS/PL considerations, and testing activities for assessing the reality and miti¬ 
gation of various technical issues associated with using PLs. 
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II. Pseudolite Signal Design Considerations 

Pseudolites operating within the GPS frequency bands (L t : 1565-1585 MHz 
or Li. 1217-1237 MHz) can be configured to serve a limited area with a power 
level low enough to preclude appreciable interference to standard GPS signals. 
One application is the PL bubble concept proposed for operating near the ends 
of airport runways to augment real-time kinematic positioning based on GPS 
carrier-phase measurements 3 (also see Chapter 15, this volume). 

For PLs designed to cover a larger area, such as an entire airport or terminal 
area, potential interference to GPS signals is a key technical issue. Although 
operation outside the GPS bands is a possibility, this option ultimately adds rf 
front end complexity and cost to user equipment. 

In this chapter, a PL signal structure is described that can operate within the 
L band and mitigate or virtually eliminate the near-far issue. First, however, a 
brief review is given of previous implementations used on DOD test ranges and 
other approaches proposed for civil applications. 

A. Previous Pseudolite Designs 

Table 1 compares two DOD implementations and three other proposed 
approaches in terms of the diversity of techniques employed for interference 
mitigation: signal (code, frequency, time) and spatial separation. Pseudolites, 
called ground transmitters (GTs) during the Phase I GPS program, were used to 
augment GPS for testing user equipment at Yuma, AZ before there were enough 
satellites for navigation. 4 The PL signal structure was the same as used on the 
GPS satellites, except for the data content. In those days, interference to the 
satellite signals was avoided by using different gold codes (Code Division Multi¬ 
ple Access, or CDMA) and keeping the user equipment under test at an adequate 
distance from the four GTs to prevent dynamic range sensitivity problems. The 
same was true in more recent Space Defense Initiative (SDI) testing with GTs 
built by Stanford Telecom (STel) to augment GPS for DOD’s Range Applications 
Program. However, a Time Division Multiple Access (TDMA) signal (on-off 
scheme) was also implemented to prevent interfering with a co-located GPS 
receiver that provided real-time synchronization to the GPS satellites. At transmit 
times, this receiver simply blanked out the signal. 

The use of PLs for civil aviation applications was first proposed in 1984. 5 
The “near-far” problem was recognized and addressed at that time, along with 
suggested solutions, but not carried any further. Three signal diversity techniques 
suggested as a solution were as follows: 

1) Pulsed signals with random or fixed cycle rates, a TDMA variation 

2) Signals transmitted at a frequency offset from L , (1575.42 MHz), but 
within the same frequency band as GPS, a variation of frequency division multiple 
access (FDMA) 

3) Alternative codes that have a longer sequence than the existing GPS codes, 
a variation of CDMA 

The first two alternatives were favored, with a preference toward the first. 
Although technically viable, the FDMA approach, including the possible use of 
L (1227.6 MHz) or frequencies outside the GPS band, was considered more 
costly given the state of GPS receiver technology at the time. The CDMA 
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Table 1 Pseudolite interference mitigation techniques—previous designs 


Pseudolite signal 
implementation (I) 
or 

recommendation 

(R) 

Interference mitigation techniques 

Code Frequency Time 

division division division Pseudolite 

multiple multiple multiple spatial 

access access access separation 

(I) Yuma ground 
transmitters for 
the Phase I GPS 
tests (1977) 

Different gold None None Unspecified 

codes 

(I) Stanford 

Telecom ground 
transmitters for 
Space Defense 
Initiative (1989) 

Different gold None Pulsed (simple Unspecified 

codes on-off) 

(R) Reference 5 

Longer U or Pulsed at random Unspecified 

sequence L, ± 15 MHz or fixed rate 

codes 

(R) Reference 6 

Different gold None Pulsed with: > 30 km 

codes 10% duty cycle, 

random pattern 

(R) Reference 7 

Different gold L\ ± A/ Pulsed with: No constraint 

codes (Afa 30 kHz) 1:11 duty cycle, 

fixed pattern, 
fixed offset 
between PLs 


approach with different pseudorandom noise (PRN) codes could be part of the 
diversity solution, but longer sequence codes would not add significant margin 
against cross-correlation interference. 

As part of the Radio Technical Commission for Maritime (RTCM) user activi¬ 
ties, a more definitive pseudolite signal structure was proposed in 1986. 6 All 
three of the above multiple access techniques were considered, but pulsed TDMA 
was the only approach recommended, because it made the least impact on the 
design of GPS receivers based on the state of technology at the time. Subsequently, 
flaws in that TDMA scheme were observed with respect to a class of “nonpartici¬ 
pating” receivers, some of which are still in use today. This led to a modified 
TDMA scheme, which was proposed in 1990. 7 Despite these proposals, fear of 
the near-far problem remained, and rightly so, because a limited interference 
margin can still exist with only code (CDMA) and pulsing (TDMA) employed. 

B. New Pseudolite Signal Design 

Fortunately, GPS receiver technology has advanced to a point that the FDMA 
approach is now viable, which improves the solution to the near-far interference 
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issue significantly. As a result, a more effective signal structure has been proposed 
that combines: good C/A codes a frequency offset that takes advantage of the 
code cross-correlation properties, and a good pulsing scheme. 8 

/. C/A Codes 

A recent study identified 19 C/A codes that were considered best (in terms 
of cross-correlation level) of the 1023 possible codes in the GPS C/A-code 
family. 9 These codes were selected for use with the FAA’s wide area augmention 
system (WAAS), which will augment GPS with a ranging signal transmitted at 
L, from geostationary orbit. Although the 19 codes selected were the best, there 
remain 712 balanced codes, almost as good, that could be used for PLs. (See 
also Chapter 3 in the companion volume for additional discussion of Gold codes 
and cross-correlation properties.) 

2. Frequency Offset 

The proposed offset is 1.023 MHz on either side of L , at 1575.42 MHz, which 
places the PL carrier in the first null of the GPS satellite C/A-code spectrum. 
This offset is much more effective and one that at least some current GPS 
receivers can accommodate, because it resembles a large Doppler frequency 
offset. Furthermore, it virtually eliminates any code cross-correlation with the 
GPS satellite signals. This spectrum is shown in Fig. 1, which is a typical spectrum 
of a GPS C/A code [PRN 2]. It was reported years ago that the cross-correlation 
between C/A codes at different frequencies was simply the magnitude of another 
C/A-code spectral line at the offset frequency. 10 This property has since been 
verified. Because all the codes have similar spectral characteristics, the spectrum 
of any C/A code tells the story. As indicated in Fig. 1 and discussed further in 
Appendix A, the spectral lines near the null are well below —80 dB and are 60 
dB lower than those at a zero frequency offset. This property has been used to 



Fig. 1 Spectrum of the GPS C/A-code PRN 2. 
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advantage in the Global Orbiting Navigation Satellite System (GLONASS) signal 
structure." All the GLONASS signals have the same maximum length codes 
with the first nulls at 511 kHz offsets from the carrier. The GLONASS FDMA 
signals are spaced at 562.5 kHz intervals, 51.5 kHz from the nulls of the adja¬ 
cent signal. 

Recently, it has been shown that cross-correlation levels can be somewhat 
higher than indicated above when the code chip boundaries are not aligned. 12 
However, when the carrier frequency is offset by 1.023 MHz, and the carrier/ 
code frequency ratio of 1540 is maintained, the code of the PL is shifting with 
respect to the satellite codes at a rate of over 664 chips per second. Thus, any 
cross-correlation is noise-like and averages to a lower rms level. As pointed out, 
this is still interference. 12 An in-band frequency offset of 1.023 MHz lowers the 
rms interference by about 8 dB, but does not eliminate it. However, it does 
eliminate cross-correlation problems. 

3. Pulsing Scheme 

The FDMA approach described above is a variation of an earlier proposal, 
which suggested a continuous signal with a longer code on the fringe of the 
allocated L\ band at 1560 or 1590 MHz. However, a continuous in-band PL 
signal could still cause interference problems, depending upon its power level 
compared to GPS satellite signals. For example, a PL signal whose power is set 
to be received at 37 km (20 n.mi.) will be 60 dB stronger at 37 m (0.02 n.mi.), 
which will be approximately 30-40 dB above the noise. Thus, although offsetting 
the frequency is a good idea, TDMA pulsing may still be required to avoid this 
situation. Furthermore, if TDMA is used, there is no reason to operate on the 
fringe of the L t allocated band, because that large a frequency offset and a longer 
code are both undesirable due to the additional burden on the GPS receiver, even 
with modem technology. 

With a good pulsing scheme, the impact on the reception of GPS signals can 
be made essentially transparent. The receiver treats it as a continuous signal, 
provided that it is designed to suppress pulse interference, as most modem 
receivers are, even if by accident. As noted in the RTCM work, 6 any hard-limiting 
receiver or “soft-limiting” receiver will clip the pulses and limit their effect, but 
still pass more than enough pulse power to track the PL signal itself. A “soft- 
limiting” receiver clips the incoming signal-plus-noise at two to three times the 
rms noise level, resulting in more, but still negligible degradation in signal-to- 
noise performance. 

All modem digital receivers are either hard-limiting or possess the soft-limiting 
property through precorrelation quantization. Although this is not true for the 
older, analog military receivers, their wideband automatic gain control (AGC) 
suppresses the pulses for the same effect. Any cross-correlation problems in 
those receivers caused by PL transmissions are eliminated with the proposed 
frequency offset. 

The pulsing scheme presented here differs from previous proposals. 67 It 
accounts for the fact that the PL message symbol rate could be 50 N sps (where 
N = 1, 2, 4, 5, 10, or 20). The recommended pulse pattern is illustrated by the 
Xs shown in Fig. 2. These repeat every 11 ms, and thus, would never be synchro- 




tu tiWt fmmhrbmpaaimdirJif Purchased from American Institute of Aeronautics and Astronautics 


56 B. D, ELROD AND A. J. VAN DIERENDONCK 


n 


C/A Cod* Cycle* 

□ 

□ 

□ 

□ 

□ 

□ 

O 

□ 

□ 

X 

X 

£££ 

£ 

£ 

£ 

£ 

£ 

£ 


£ 

H 

1 

u 




a 



z 




a 

Z 



U 


— 

“ 


~ 


[X 


2 


li 




u 


L 




■ 

u 

■ 



u 

■ 

■ 




■ 

Li 

3 


■ 

u 




a 







LI 




u 



— 




4 




u 



■ 

U 







LI 




u 

■ 

z 




5 





u 




U 







Li 




a 

■ 




6 






LI 




V 







u 




Y 




7 







u 




u 







LI 




U 

■ 


8 

u 






1 

u 




u 







LI 



u 

U 


9 


a 







u 




u 







u 


_ 


Y 

10 


■ 


■ 






u 



■ 

U 







Li 




11 


r 

mu 







u 




vl 







LI 



X = PL, Pulses, Y = PL, Pulses, etc. 


Fig. 2 Illustration of pseudoiite pulse pattern. 


nous with a received GPS bit pattern. Each code cycle (1 ms) is divided into 11 
slots, each with a width of 1/11,000 s (90.90909 p,s). The pulses would never 
be synchronous with a received GPS bit pattern, as would occur with the 
RTCM design. 6 

Only one of these slots contains a pulse, so the duty cycle is 1/11. There will 
be 20 IN pulses per symbol, but only one is required because the receiver would 
integrate the energy over the entire symbol period. Because every slot is filled 
once every 11 ms, the entire C/A code would be received every 11 ms. The 
clocking rate for the slots is 1/93 of the C/A-code chipping rate of 1.023 MHz, 
or 11 kHz. It is noted in Appendix B that with a duty cycle of 1/11, the loss 
of GPS satellite signal-to-noise ratio, in either an analog receiver with pulse 
suppression, or a digital receiver with natural soft-limiting, is less than 1.5 dB 
when close to the PL. 

Potential mutual interference when multiple PLs are installed in an area was 
reported in the RTCM work. 6 With the RTCM pulse pattern, a minimum distance 
between PLs would be required, because the transmission of pulses simultane¬ 
ously from each PL could result in the simultaneous reception of multiple pulses. 
But, that need not be, because the pulse timing of multiple PLs can be offset. 
Unfortunately, the RTCM pattern had irregular times between pulses, so that 
there would still be some simultaneous receptions. 

However, with the pulse pattern shown in Fig. 2, a suitable pulse-timing 
offset would prevent simultaneous reception from ever occurring, except at large 
distances where the PL signal would be of little consequence. For example, 
consider two PLs where the pulse timing was offset by 4 ms as indicated by the 
Xs and Ys in Fig. 2. Because the minimum transmit time separation between X 
and Y pulses is 4/11 ms, only receivers with a differential range more than 110 
km from the PLs would encounter simultaneous reception. At that distance, the 
received powers of either PL would be negligible. 
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4. Pseudolite Carrier Tracking 

A potential misconception about pulsing a PL signal is that it will prevent 
the tracking of PL carrier phase. This is not true if the pulses occur at a high 
enough rate, which is the case in the proposed scheme. One or more pulses will 
always be integrated along with noise over each symbol. The result is transparent 
to the tracking loops, because the phase change due to Doppler uncertainty over 
a symbol is negligible. To the tracking loops, it looks like a continuous signal. 


5. Pseudolite Transmit Power 


Given an average received power P r through a receiving antenna (with gain 
G a ) at a distance d (in n.mi.) from a PL, the average transmitted power (P,) at 
Li ± A/ is as follows: 


P r + 20 log, 0 -- G a 




— P r + 20 log io d + 101.75 — G a 

where X., = 0.00010275 n.mi. is the signal wavelength corresponding to L,. As 
an example, consider an average received power of —130 dBm at 20 n.mi. (37 
km) through an antenna gain of —10 dB (assumed for a small negative elevation 
angle). The average transmitted power is 7.77 dBm, or about 6 mW. The peak 
power for a duty cycle of 1/11 is then 18.18 dBm, or about 66 mW. 


III. Integrated Differential GPS/Pseudolite Considerations 

The introduction of PLs has two key objectives: signal augmentation and data 
link enhancement. The first is to increase the number of available signals and, 
thereby, improve or maintain the geometrical quality for user position determina¬ 
tion. The second is to provide an integrated capacity for supplying key data to 
users for GPS (and PL) integrity warning and differential corrections to improve 
positioning accuracy via code-based local differential GPS (LDGPS) and poten¬ 
tially carrier phase-based kinematic differential GPS (KDGPS) techniques. 

Although the definition of a compatible signal format is one of the most 
critical requirements for meeting these objectives, there are other implementation- 
related aspects that need consideration as well. This section addresses five of 
these: PL siting, PL time synchronization, antenna location on a user aircraft, 
the PL data message, and the filter algorithm for integrated GPS/PL measure¬ 
ment processing. 


A. Pseudolite Siting 

It is well known that the GPS geometry (quantified in terms of HDOP and 
VDOP factors) will vary over time and user location, even with a full 24-satellite 
constellation 13 (also see Chapter 5, the companion volume). It is also known that 
at times significant degradation (VDOP >> 6) can occur if fewer satellites are 
active. The utility of PLs for geometric enhancement lies in the fact that lower 
DOP values with less temporal variation can be achieved. 
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Figures 3 and 4 illustrate representative HDOP and VDOP profiles determined 
for a full constellation 13 and with one satellite inactive for cases where the user 
employs the best four GPS satellites or “all-in-view” with augmentation by one 
or two PLs. 2 In this illustration, the user was assumed to be at an altitude of 200 
ft (the decision height for a Cat. I approach) with the PLs located 1 n.mi. ahead 
and/or 2 n.mi. behind. It is evident that in this situation, a two-PL augmentation 
would enable HDOP ^1.0 and VDOP ^ 1.5 continuously. At those levels, code¬ 
based LDGPS with the capability to correct uncorrelated pseudorange errors to 
< lm (2cr) would meet the current Cat. Ill (horizontal) and Cat. II (vertical) 
sensor accuracy requirements of 4.1m and 1.4m (2 ct), respectively for aircraft 
precision approaches. 14 With the addition of ranging-capable satellites in geosta¬ 
tionary orbit (e.g., via the FAA’s WAAS implementation) to offset one or two 
inactive GPS satellites an even lower VDOP is possible with favorable implica¬ 
tions for meeting the most stringent Cat. Ill (vertical) accuracy requirement of 
0.4m (2 ct). 14 

An alternate approach to meeting the Cat. Ill vertical accuracy requirement 
for precision approaches is based on KDGPS techniques with PL augmentation. 
Recent flight tests have demonstrated that a varying user/PL geometry will 
permit rapid resolution of carrier-range ambiguities to provide KDGPS-based 
measurements with centimeter-level precision 3 (also see Chapter 15, this volume). 
The configuration adopted for this technique utilizes a bubble concept in which 
two PLs located ahead of the runway threshold on each side of the glideslope 
transmit a standard, unpulsed L r CI A signal at a low power level sufficient to be 
received within a hemispherical bubble centered at each PL. A descending aircraft 
would enter the common bubble zone, acquire the PL signals, resolve GPS carrier- 
phase ambiguities, and emerge to continue the approach with KDGPS-based 
corrections supplied via a separate (non-PL) datalink. As proposed, four active 
PLs per runway would be needed to support landings in either direction. 

On the other hand, wider coverage large bubble PLs with interference mitiga¬ 
tion inherent in the signal design could provide service to multiple runways over 
an entire airport region. Location of the PL antenna off runway and at suitable 
elevation would also be desirable from the standpoint of user antenna requirements 
(discussed in Sec. II.C). Pseudolite-siting requirements (number and location) to 
accommodate multirunway situations, user antenna considerations, and local 
constraints are key issues needing further analysis and testing. 

B. Pseudolite Time Synchronization 

Synchronizing a PL’s clock to GPS time can be achieved in two ways—one 
in which a PL is collocated with an LDGPS reference receiver (RR), and one 
in which it is remote from the RR that is tracking its transmitted signal. In the 
latter case, the RR sends corrections to the remote “slave” PL to correct itself 
and also sends message data (e.g., code- and/or carrier-based DGPS corrections) 
to be transmitted. These two different PL configurations are illustrated in Fig. 
5. In the co-located configuration, the RR shares the transmit/receive antenna 
with the PL, which also allows for self-calibration. 

The type of configuration used would depend upon whether or not there is 
more than one PL at a given local region. If there is only one, the collocated 
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Fig. 3 Worst-case dilution-of-precision profiles at 35°N with and without PLs for 
a full (24) GPS constellation. 
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Fig. 5 Master and slave pseudolite configurations. 


approach is more desirable, especially if line-of-sight (LOS) visibility to the RR 
might be a problem. If there is more than one PL, the remote approach with a 
common RR for synchronization is more desirable. However, if LOS visibility 
problems exist, having a RR collocated with each PL would enable time synchro¬ 
nization via GPS common-view, time-transfer techniques. However, only one 
RR can be used for deriving satellite LDGPS corrections. In this case, master 
and slave PLs would be designated with the master clock and LDGPS corrections 
derived centrally and distributed to each remote PL, which would slave its time 
to the master RR clock. This approach is similar to that used by various DOD 
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test and training ranges coordinated by the Range Applications Joint Program 
Office (RAJPO), for which P-code PLs were developed by Stanford Telecom. 

1. Master Pseudolite Configuration 

Both the RR and the PL signal generator derive their timing coherently from 
the same stable frequency standard. The signal generator/pulser electronics mod¬ 
ule pulses the transmission of the PL signal in order to minimize interference to 
both participating and nonparticipating GPS receivers, as discussed in Sec. IB3. 
This multiplexing also allows the RR to receive the satellite signals via the same 
antenna. In fact, by providing a suitable calibration path, the RR can also track 
the output of the signal generator. In this way, the collocated PL is self-calibrating, 
and the transmitted PL signal will be synchronized to the same clock that is used 
to derive the differential GPS corrections. This is true even in the multiple-PL 
scenario, where the slave PLs receive differential corrections from the master 
PL. In this case, the time solutions of the slave PL receivers will be referenced 
to the master PL’s clock, because the differential corrections are computed with 
respect to that clock. 

2. Slave Pseudolite Configurations 

The slave PLs need not have receivers if they can be tracked by the RR. The 
RR simply supplies corrections to the PL for correcting its clock via a number- 
controlled oscillator (NCO) and provides the satellite differential corrections for 
modulation on the PL signal. Because the RR can update the PL continuously, 
a good quality crystal oscillator will suffice as its frequency reference. Otherwise, 
its configuration is identical to the co-located configuration but without the RR 
and the self-calibration path. 

C. User Aircraft Antenna Location 

Reception of the PL signal by a user will be affected by the aircraft antenna 
location and the PL position relative to its approach path. Ideally, the PL signals 
would be received by a top-mounted antenna, the same one used for receiving 
the satellite signals. This could be accommodated by locating the PL antenna at 
an appropriate elevation and offset from the glideslope. Even if the line of sight 
to the PL is below the aircraft antenna horizon, the increased signal level in the 
vicinity of the PL will tend to cancel the loss in antenna gain. 

If the aircraft were to pass directly over a PL, a larger angular gradient would 
be available to support the positioning process. However, a bottom-mounted 
antenna and a separate front end to interface the antenna(s) to the receiver may 
be needed. A configuration with a bottom-mounted antenna would likely have 
more sensitivity to ground-reflected multipath and interference. A top-mounted 
antenna is less sensitive, because reflective surfaces on the aircraft are typically 
small relative to the C/A-code chip width (293 m). In addition, the use of a 
bottom-mounted antenna would introduce additional sensitivity to lever arm 
uncertainty and knowledge of aircraft attitude information. On the other hand, 
it could be used directly, if a PL were to be employed only as a datalink. A 
bottom antenna is required for the small bubble concept. 3 Obviously, further 
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analysis and testing are needed before the aircraft antenna location issue is 
resolved, one way or the other. 

D. Pseudolite Signal Data Message 

A PL offers the possibility for an order-of-magnitude increase in the data rate 
(up to 1000 bps vs 50 bps for GPS) that can be supported via a GPS-compatible 
signal with essentially a firmware change in the user receiver. Validation of data 
link performance at this rate is a key test objective. A closely related issue is 
the PL message capacity required to support GPS and PL integrity updates and 
DGPS corrections (code and/or carrier). 

Some have proposed that for the latter, a 2400 bps data link would be needed 
to supply raw carrier-phase and other GPS information from a reference receiver 
at a 0.5-1.0 Hz update rate. Efforts by RTCM SC-104 (Special Committe 104) 
and others contend that KDGPS could be supported with much lower data transfer 
requirements (< 1 kbps) based on the use of carrier phase corrections, not raw 
data. 15 A related issue for testing is whether PL aiding of carrier-phase ambiguity 
resolution and cycle slip maintenance procedures could reduce the data require¬ 
ments still further. 

1. General Format 

For efficiency, the general format is patterned after the WAAS format 16 with 
three differences—it may or may not include forward error correction (FEC), a 
7-bit distributed time word is added and the data rate can be higher than 250 
bps, up to 1 kbps without FEC. This general PL format is shown in Fig. 6. The 
time word is added as a convenience to the user receiver, because PL time is 
not the same as GPS time. This is because the code-chipping rate is offset in 
frequency by approximately ± 664 chips per second, the feature that eliminates 
cross-correlation with the GPS signals. The actual PL frequency can be chosen 
(at an offset from the spectral null of 710.4166667 Hz) so that the PL week is 
exactly 393 s shorter or longer than the GPS week. In fact, PL time is, at any 
GPS time, as follows: 

393 

i ± _ _ 

604,800 

depending upon which null is selected. The PL time is distributed over three 
250-bit subframes making up a 21-bit word. The 21-bit time word represents the 





4 -Direction of Data Flow From PL; Most Significant Bit (MSB) Transmitted First 



24-Bits 

, | | 205-Bit Oata FiaM 1 p arity 




-7-Bit Distributed Time Word 

-6-Bit Message Type Identifier (0-63) 

-8-Bit Preamble of 24 Bits Total in 3 Contiguous Blocks 


Fig. 6 Pseudolite message format. 
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PL time of week at the. start of the currently transmitted 24-bit preamble, starting 
with 0 at the beginning of the week. This PL time word also serves as the 
reference time for the data in the messages. 

The 24-bit parity is the same as the cyclic redundancy check (CRC) parity 
specified for the WAAS. 16 The data field consists of 205 information-bearing bits. 

2. Message Types 

To avoid confusing the message types with those of the WAAS, the message 
type numbers start at 40. Table 2 lists a tentative set of PL message types. Every 
message would include a certain number of GPS/PL signal integrity flags to 
provide a short time to alarm capability. 

3. Message Content 

The contents of the messages are somewhat different from the message contents 
defined for the RTCA SCAT-I (Special Category I) DGPS. 17 This is to accommo¬ 
date the 205-bit data fields and to provide data that are more consistent with 
Category II and III precision approach requirements. The SCAT-I messages are 
quite inefficient in requiring too much signal bandwidth. Preliminary message 
contents show more than adequate bandwidth in using the maximum 1 kbps 
capability of the proposed PL signal. Every message broadcast contains integrity 
flags for several PRN numbers, with a minimum of 11, including the transmitting 
PL. This allows for a positive integrity indication at least once per 0.5 s. PRN 
numbers that accompany the flags do not have to be broadcast in order, so an 
alarm can always be inserted in any 0.25 s message providing a maximum time 
to alarm of 0.5 s. 

Message Type 42 consists of a slight modification to the content of the 
corresponding RTCM carrier-phase corrections message. 15 Each can accommo¬ 
date four satellites/pseudolites. Therefore, both pseudorange and carrier-phase 
corrections can be broadcast. 

E. GPS/Pseudolite Navigation Filter Considerations 

When pseudorange (or carrier-range) measurements are processed by the user 
navigation filter, modeling to accommodate the nonlinearity in the measurement 


Table 2 Pseudolite message types 


Type 

Contents 

40 

Don’t use this PL for anything (PL testing) 

41 

Integrity flags/Pseudorange corrections 

42 

Integrity flags/Carrier-phase corrections (if 
required) 

43 

Integrity flags/PL location and PRN assignment 

44 

Integrity flags/PL almanacs 

45 

Integrity flags/Precision approach path 
definition 

46 

Integrity flags/Special message 

63 

Null message—Alternating Is and 0s 
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model can be ignored in the case of satellites, but not necessarily for PLs. As 
the user range to the PL diminishes during an approach, the impact of the 
nonlinearity is to introduce an apparent bias into the measurements. If this error 
is comparable to the measurement error, then a standard extended Kalman filter 
(EKF) will yield inferior performance. Filter divergence may occur, as the combi¬ 
nation of measurement and nonlinearity error exceeds the filter’s own computation 
of rms error, and it rejects new measurement data. 

One approach to preventing filter divergence of this sort is to choose suffi¬ 
ciently large a priori measurement variances that include worst-case nonlinearity 
effects. Unfortunately, this requires identification of what is worst case. More 
importantly, high constant measurement variances will cause sluggish perfor¬ 
mance at other times when there is a negligible nonlinear effect present. 

Another approach, as outlined in Appendix C, is to increase the filter sophistica¬ 
tion by using a Gaussian second-order (GSO) filter that is similar to a (linearized) 
EKF but includes a quadratic component of the measurement nonlinearity. The 
key benefit is that it offers improved performance when the nonlinear effect is 
present but reverts naturally to a standard EKF when it is not. A possible downside 
is that more software and processing time are required, although some approxima¬ 
tions are possible depending upon the scenario. Pseudolite siting relative to the 
approach trajectory will also be a factor in this. 


IV. Pseudolite Testing 

The PL concept for augmenting GPS, as described above, is intended to provide 
users with the following potential benefits: 1) enhanced local area navigation 
performance via integrated GPS/PL positioning using single or multiple PL 
signals with airportwide or even terminal area coverage depending upon the 
power level and antenna configuration (s) employed; 2) more data link capacity 
to support DGPS operations directly (code and/or carrier-based) at a multiple N 
of 50 bps (where N = 2,4,5, 10, or 20); and 3) mitigation of potential interference 
to standard GPS signals through the code, frequency, and time diversity techniques 
employed in the PL signal design. 

Initial testing to verify these features has been conducted under a research 
and development project sponsored by the FAA, 18 and further comprehensive 
testing is planned. The following subsections describe the test results. 

A. Pseudolite Interference Testing 

Initial laboratory tests to evaluate PL interference impacts on GPS receivers 
was conducted by applying simulated GPS and PL signals at various power levels 
to several commercial receivers. After initial signal acquisition and steady-state 
operation by the GPS receiver was reached, a PL signal was injected at gradually 
increasing power levels. Signal quality reported in terms of GPS carrier power- 
to-noise power density ratio ( C/N 0 ) was recorded as a function of average PL 
signal power. The plots in Fig. 7 show the test results obtained for two receivers 
under the following PL signal conditions: 1) no pulsing, no frequency offset; 2) 
no pulsing, frequency offset applied; 3) pulsing applied, no frequency offset; and 
4) pulsing and frequency offset applied. 
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Fig. 7 GPS degradation vs pseudolite signal power with and without pulsing and 
frequency offset. 


With no pulsing or frequency offset applied, the degradation is significant at 
just a 20 dB advantage in PL over GPS signal power. With either pulsing or 
frequency offset applied, the degradation is less, although pulsing is the more 
effective feature. With both applied, the degradation in reported C/N 0 from the 
no PL signal case is small (s 2 dB) over a range of 60 dB or more in the ratio 
of PL to GPS signal power. The results for this case from testing four receivers 
(two L r C/A only and two LJL 2 cross-correlating types) are shown in expanded 
scale in Fig. 8. When compared to a user receiving a PL signal at a range of 20 
n.mi. (37 km) at a level comparable to GPS, this means that a degradation of 
only 2 dB would be experienced by a user only 0.02 n.mi. (37 m) away but 
receiving the PL signal at 60 dB higher power! 

Note that receiver A’s degradation with a continuous signal, or with either a 
frequency offset or pulsing, but not both, is more than that of receiver B. This 
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Fig. 8 GPS degradation vs pseudolite signal: 4 receivers with pulsing and 
frequency offset. 
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illustrates the difference between soft-limiting and hard-limiting. Receiver A 
uses multibit sampling. However, with both features turned on, receiver A’s 
performance is quite good. Note also that receiver B’s performance tails off as 
the PL signal gets quite strong (> 50 dB above GPS). Receiver B is a discontinued 
model, and receiver C is its replacement. Receiver C did not exhibit the same 
behavior, which was probably because of slow saturation recovery in the receiver 
B’s front end. For receivers with fast enough saturation recovery, the results 
agree well with Eq. (B12) and Eq. (B13) of Appendix B. 

The PL interference test results discussed above were obtained with high-end 
receivers capable of reporting signal quality data. Testing of other relatively low- 
cost commercial receivers was also done on a qualitative basis. These receivers 
were placed at varying distances from an antenna radiating a PL signal with 
pulsing and frequency offset features applied or not. With the PL signal off, each 
receiver was set up for normal GPS operations. With the PL signal on at peak 
power, each receiver was moved toward the PL antenna, and the separation 
distance was recorded when anything anomalous appeared on its display. The 
results for the four units tested in this manner showed no effect until within 1-5 
m of the PL antenna at a peak radiated power of +15 dBm. 

B. Pseudolite Data Link Testing 

Tests were performed to verify the capability to transmit, receive, and demodu¬ 
late the message data encoded on a PL signal with the pulsing, frequency offset, 
and higher data rate features described above. During laboratory tests at Stanford 
Telecom, a signal generator/pulser unit was interfaced directly to a GPS/PL- 
capable receiver. Pseudolite signals with a fixed data message were transmitted 
at different power levels to simulate operations at various PL/user ranges. Results 
of these initial tests indicated virtually error-free data reception, with the received 
signal quality ( C/N a ) at a level of 35 dB-Hz or more. 

To support planning for more comprehensive PL testing and evaluation, prelim¬ 
inary flight trials of a DGPS data link provided by one PL were arranged. A 
GPS/PL reference receiver and the signal generator/pulser unit used for the 
laboratory tests were installed at the FAA Technical Center (FAATC) and config¬ 
ured to transmit the PL signal from an antenna on the hanger roof. The test 
aircraft was an FAA Aerocommander (N50) with a top-mounted, low-profile 
GPS antenna, a GPS/PL-capable receiver, and data processing/storage equipment. 
For test purposes, the PL message data (LDGPS) corrections) were encoded in 
a 250-bit WAAS message format and transmitted without forward error correction 
at 250 bps. 

Flight profiles flown with the test aircraft included straight-in approaches to 
Runway 6 at FAATC and various holding patterns ranging up to 29 km (15 n.mi.) 
from the PL. Received PL message data and signal quality data (C/Nd) were 
recorded onboard during flight segments, such as A, B, and C in Fig. 9. 

Corresponding CIN 0 profiles and occurrences of individual data message errors 
are shown in Fig. 10. The PL data error occurrences appear to be predominantly 
associated with PL antenna gain reduction/obstruction phenomena during turns. 
Future flight trials will help address the impact of aircraft antenna type on PL 
signal-tracking performance. The test aircraft will be equipped to receive PL 
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Ground track 




Fig. 9 Examples of flight test segments at the FAA Technical Center during pseu- 
doiite data link tests. 


signals from top- and bottom-mounted, low-profile antennas, and a high-profile 
(blade) antenna. These trials will also assess PL data link reliability at other data 
rates (500 and 1000 bps). 


C. Navigation Performance Testing 

The initial laboratory tests and preliminary flight testing have focused on 
demonstrating the feasibility of PL tracking and interference mitigation to GPS 
signals with the pulsing and L| offset features applied. Additional aspects that 
remain to be assessed include the following: 

1) PL code and carrier measurement quality 

2) accuracy performance enhancement in LDGPS/PL and KDGPS/PL modes 

3) sensitivities to PL siting (number/location) and user antenna (type/location) 

Future flight trials are planned to assess navigation performance in LDGPS/ 

PL and KDGPS/PL modes, and the sensitivity to PL siting and type/location of 
user antenna employed. 

Independent KDGPS-based tracking of the test aircraft will be used as a 
truth source for assessing the consistency of PL ranging data (code and carrier) 
throughout a terminal area flight envelope. 




C/No (dB-Hz) C/No (dB-Hz) C/No (dB-Hz) 
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Fig. 10 Pseudolite signal profiles during data link tests. 
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Appendix A: Interference Caused by Cross Correlation Between 

C/A Codes 


The effect that one C/A code has on another with respect to cross-correlation 
properties was described in Ref. 10 (see also Chapter 3, the companion volume). A 
more rigorous derivation is provided here based on the receiver signal-processing 
model shown in Fig. Al. 

In this signal-processing model, the input signals at G>, in terms of in-phase 
and quadraphase components, are the following desired signal: 

/,i(0 = A,C, (f)cos(2n A/r + <(.,) 

Quit) = A,Ci(r)sin(2irA/;r + (fr) (Al) 

and the undesired signal 

Iji(t) = AjCj{t 4- t ; (0)cos[2tt(8/ + Afi)t + <]>,] 

= AjCj(t + T,(r))sin[2ir(S/J + A f,)t + <j> y ] (A2) 

where 

A, , Aj = signal amplitudes 
C,(r), Cj(t) = signal C/A codes 
A/, Afi = signal Dopplers 
<t>i, 4>y = signal phases 

8 fj = frequency offset of undesired signal 
Tj(t) — time offset between signals 

The Doppler removal process eliminates the desired signal’s Doppler and phase 
so that the signal components at ® in Fig. Al are as follows: 

I a (t) = A,Ci(t) (A3) 

Qdt) = 0 


IjiU) = AjCjit + T ; (f))cos[2Tr(8/y + A fi - A f)t + 4> y - 4>J 
Q fl ( 0 = AjCjit + T y (r))sin[2rr(8^ + Afi - A fi)t + <h, - 4>,] 



Fig. Al Receiver signal-processing model. 
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where the Doppler difference is as follows: 

A f„ — bfi + Afi - Af (A5) 

For simplicity, assume that the code transitions line up. Thus, for the moment, 
assume the following: 

TjU) = NT C (A6) 

where N is an integer, and T c is a chip width (1/1,023,000 s). Also, assume full 
correlation for signal i. Then, at ® we have the following: 


/„(0 = A, 
fi 3 (0 = 0 


(A7) 


7, 3 (r) = A.jCj{t + NT :: )C;(t)cos(2iTAf :j t + 

0,3(0 = AjCj(t + NT c ')C,(t)s'm( 2rr Af ;j t + <J> ; . - <j> ; ) 


(AS) 


Now, because of the cycle-and-add property of the C/A-codes, Eq. (A8) becomes 
the following: 


/j3(0 = Aj C k (0cos(2-irA/ y / + <t>, - 4>,) 
2,3(0 = Aj C k (/)sin(2irA/^ t + <t>, - 4> f ) 


(A9) 


where C k (t) is another code in the same family, which is different for each value 
of N. The signal components at ® are then as follows: 


7,4(0 = A t T 
2,4(0 = 0 

/,4(0 = Aj [ C ft (Ocos(2'7rA/i ; -f + ^ - 4),) dr 

j 0 

2 , 4(0 = Ay f C t (0sin(2TrA/yr + <})y - 4>i) dr 

J 0 

where T is a multiple M of 1 ms C/A-code repetition periods. 
Power in the two correlated signals is given as follows: 


(A10) 

(All) 


2Pm = 7|t + Qi = AfT 2 


2Pj (Ay ) = I]< + Qj 4 


~ Aj 


\ C t (0cos(2ir/yr) dr + f Cd0sin(2rr/yr) dr 

•7 0 -7 0 


(A 12) 


(A13) 


Note that, through expansion using trigonometric identities, the dependence upon 
the phase difference has been removed in Eq. (A13), which resembles the expres¬ 
sion for a Fourier power spectrum component at the Doppler difference. The 
code repeats at a 1-kHz rate, and the integration over each 1 ms code period is 
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identical for Doppler differences of multiples of 1 kHz. Thus, the ratio of Eq. (A13) 
to Eq. (A 12) the cross-correlation power ratio, can be reduced to the following: 


w = % 


_1_ f 2L 

2L Jo 


Ct cosi 


/ITT 


t dt 


+ 1 Yl C c * si "(f') * 


(A14) 


where 


A fa = nI2L = 1000 n Hz 
2 L = 0.001 s = TIM 


The term in the brackets of Eq. (A14) can be recognized as the nth power 
spectrum component. Obviously, Eq. (A 13) can take on values for Doppler 
differences other than multiples of 1 kHz, depending upon the value of T, but 
peaks at the 1 kHz lines. That is, there is cross-correlation at other Doppler 
differences, but less than at the 1 kHz lines. Because of the navigation message 
databits, M is limited to 20. If T were allowed to be infinite (very large), the 
cross correlation at these other values would approach zero. For example, if T 
was a large multiple of Doppler difference cycle periods, integration over each 
cycle would pick up a different portion of the code, or multiple codes plus a 
different portion of the code. If we take a Doppler difference of 50 Hz and 
integrate over 20 ms, then integration over the last 10 ms would cancel that over 
the first 10 ms, because the codes would simply be flipped. This is true for any 
Doppler difference that is a submultiple of 1 kHz and an integration time that 
is an integer number of Doppler difference cycle periods. For other Doppler 
differences, other than the multiples of 1 kHz, a longer integration period is 
required for eventual cancellation. 

Equation (A 14) was evaluated for equal amplitude signals for all relative code 
phases as defined in Eq. (A6) for two specific GPS C/A-code pairs. The first 
pair (PRN6/PRN28) is considered to be the worst pair of GPS codes; whereas 
the second pair (PRN7/PRN201) is made up of the best GPS code and the best 
selected for Inmarsat-3. 9 For each pair, the maximum spectrum components over 
all possible integer code phases (0 < N < 1022), are plotted in Fig. A2 for 
Doppler differences up to 20 kHz. The maximum spectrum components are also 
plotted in Fig. A3, for Doppler differences in the range 1017-1029 kHz, which 
corresponds to one of the signals transmitting in the first null of the other. This 
± 6 kHz Doppler range represents the maximum expected from satellite and 
user motion. The average over all N computed for the PRN6/PRN28 pair is also 
plotted in Fig. A3 for comparison. The average is approximately equal to the 
following spectral line envelope: 


S(f) 


1 sin z (mr/1023) 

1023 (mr/1023) 2 


(A 15) 


which is about 8-9 dB below the maximum values. Note that in both figures, it 
doesn’t seem to matter which code pair is used when determining the worst-case 
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Fig. A2 Maximum cross-correlation spectral components for two code pairs in 
Doppler difference range of 0-20 kHz. 



Fig. A3 Maximum cross-correlation spectral components for two code pairs in 
Doppler difference range of 1017-1029 kHz. 


magnitudes. This is because each of the 1023 code phases (0 s N < 1022) 
results in a different code. Consequently, codes with bad (i.e., large) spectral 
components are generated in each case, and some of these codes are unbalanced, 
as well. 

Although some of the worst-case components shown in Fig. A3 slightly exceed 
the predicted —80 dB level stated in Sec. IIB1, they are still typically below 
—70 dB. Within ±2 kHz of the null, they are, in fact, below —80 dB. More 
significantly, all components are a good 50-60 dB lower than the level near 
center frequency. This is in addition to the margin realized from pulsing with a 
1/11 duty cycle. (See Appendix B) 

The derivation presented above applies when the cross-correlation code transi¬ 
tions are lined-up, which, of course, they rarely will be. Spilker 10 points out that. 
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at the n kHz carrier frequency differences, there also exists a code frequency 
shift, which is less by a factor of 1540. Thus, for example, at the 1 kHz carrier 
frequency difference, there is a code frequency difference of 1000/1540 <** 0.65 
chips per second. Thus, the code transitions will not stay lined-up. This is 
especially true for a large frequency difference of 1.023 MHz, in which case the 
code frequency difference is approximately 664 chips per second. McGraw 12 
pointed out that cross-correlation levels can be even higher when the code transi¬ 
tions are not lined up. However, because of the rapidly changing time relationships 
between the codes, the resulting cross-correlation becomes noise-like, and simply 
becomes a noise interference. The effect of this interference is addressed in 
Appendix B. 

Appendix B: Interference Caused by Pseudolite Signal Level 

One C/A-code signal can interfere with another if it is strong enough, indepen¬ 
dent of the cross-correlation. This can certainly happen in the case of a PL signal, 
which may interfere with satellite signals as well as with other PL signals. In 
other words, a PL signal is a noise source that may jam the other signals unless 
measures are taken to mitigate this effect. The following provides an assessment 
of the impact of that jamming and a method for minimizing it; namely, by pulsing 
the PL signal. 

Background 

In general, the spreading process in a receiver’s correlator is defined as a 
signal (noise or otherwise) being passed through a filter described with a frequency 
response equal to the spectral density of the PRN code. That is, the interference 
noise density at the output of the correlator is as follows: 

No, = f S c (f)Stf) df (Bl) 

J —00 

where S c (f) is the spectral density of the reference PRN code and S/f) is the 
density of the interference or noise. The reference C/A-code has a discrete spectral 
density that can be described as follows: 

oo 

Seif) = £ C " S U" - 1000 (B2) 

where the c„ are spectral line coefficients: 8(/) is the dirac delta function; and 
the c„ vary about the envelope of Eq. (A15) in Appendix A. The reference C/ 
A-code spectral density has the property that 

S c (f) d/ = 1 (B3) 

First consider bandlimited thermal noise with density N a with a two-sided 
intermediate frequency (IF) bandwidth of B, (brick wall filter). 

N m = N 0 f' 12 Sff) d/< N 0 (B4) 

1 - s,n 

The variation of the C/A-code spectral lines averages out over the wide bandwidth. 
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A similar equation applies for wide- and narrow-band interference with a spectral 
density as follows: 

SiCfl = (B5) 

Ju Jl 

for upper and lower frequency limits /„ and/, (converted to IF) and total interfer¬ 
ence power Pi (relative to the signal power S). 

p /*min(a|,a2) 

Nat = 7~77 S,(f) df (B6) 

Ju Jl J max(-aj,a3> 

where 


a, = B,l 2, a .2 — fu ~ f/F and a 3 = /, - f, F . 

For narrow bandwidth noise interference centered at the GPS frequencies (i.e., 
somewhat less than the code chipping rate 1 IT C ), this equation becomes as follows: 

- P,T C (B7) 

This is only true for the C/A code where interference bandwidths are on the 
order of 100 kHz or greater because of the variations in the line spectrum of 
the codes. 


Pseudolite Interference 


Now assume that the interference is another C/A-code signal with the following 
spectral density: 


Srdf) = Pn. 1000 n ± 1.023 X 10 6 ) (B8) 

«—n b 


where P rL is the received PL signal power, and N B indicates the band-limiting 
effect. Because the summations in Eqs. (B2) and (B8) are over a large number 
of spectral components, it suffices to use the average envelope of Eq. (A 15), 
divided by 1000 to spread the components into a continuous spectral density, for 
evaluation. Then, Eq. (Bl) can be approximated with the following integral: 


AW 


p T 2 f 8 ' n sin 2 (TT//)sin 2 [TTC f± 1/7373 
VLc Fb ,,2 (TTfr c nir(J± \/T c )T c } 2 1 


(B9) 


This can be compared to the interference noise density from one GPS satellite 
j to another satellite i at the normal frequency given by the following: 



B i 12 sin 4 (-n//.) 
- B,n ( tt //.) 4 


df 


(B10) 


In general, Eqs. (B9) and (B10) must be evaluated numerically. This was done 
over a wide bandwidth of ±10/7/ resulting in the following relationship: 



Nqpl , , n , 

— j = 10 log,,). 



— 8.19 dB 


(Bll) 
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which indicates that the PL interferes by 8.19 dB less by transmitting in the null, 
than it would if it were transmitting at the same frequency as the satellites. 
However, this is not enough! 

The first term of Eq. (B11) can become significant as the user receiver comes 
closer to a PL. This is a key reason for PL pulse modulation, which has the 
effect of only interfering a percentage of the time equal to the pulse duty cycle, 
provided that the receiver clips the pulses. This results in a loss in received 
satellite C/N 0 of either: 


L = 


■ 10 log,, 


1 - DC + 


2 R. 


A/ 



DC 

(1 - DC) 2 



(B12) 


or 


L = 10 log, 


1 4. PpL 
3 N 0 c 


DC 


(B13) 


depending upon whether or not the PL signal is saturating the receiver’s analog- 
to-digital (A/D) sampler ( soft or hard-limiting), where Ry is the reduction in 
interference realized using a frequency offset (0.1517, if ±1.023 MHz, and 1, 
if on frequency), L m „ is the maximum sampler threshold level, L N is the one 
sigma noise level in terms of the threshold level and P pl IN q is the average received 
PL signal-to-noise density. Eq. (B12) represents the loss when the PL pulses are 
saturating the A/D and Eq. (B13) represents the loss when they are not. Note 
that in the former case, the loss is proportional to the ratio of IF bandwidth to 
the code bandwith. 

Cross-correlation can still occur, even if the pulses are clipped. However, as 
was shown in Appendix A, the cross-correlation is reduced substantially by PL 
transmission in the first null of the satellite C/A-code/spectrum (L, ± 1.023 
MHz). Pseudolite transmission at higher nulls is also possible, but this begins to 
add more complexity to a receiver designed to process both GPS and PL signals. 
Another key reason for pulsing is to prevent capturing the front end of the user’s 
receiver. This is especially important for hard-limiting receivers (1-bit samplers) 
that normally do not employ front-end AGC circuits. 


Appendix C: Navigation Filter Modeling with Pseudolite Measurements 

A GPS pseudorange measurement (and carrier range) is generally a nonlinear 
function of the satellite position vector for (r s ) and the user position vector (r„). 
This measurement is modeled by the following: 

Z = g + b u + v (Cl) 


where 


g = lr s - rj 


(C2) 



tu Mbrifr ftmfarjbmpotviWirAp Purchased from American Institute of Aeronautics and Astronautics 


PSEUDOLITES 77 

is the geometric range; b u is the user clock offset from GPS time, and v represents 
the composite of uncorrected measurement errors caused by atmospherics, satel¬ 
lite timing offsets, multipath, and noise. Given the satellite ephemeris and an a 
priori estimate of the user location (r u ), then g can be represented by the Taylor 
series expansion: 

g = go + hAr u + ArlJArJ2 + ■■■ (C3) 

where 

Ar„ = r„ — r~ (C4) 

g 0 =\r s -r-\ (C5) 

h = 3g/5r„l ru = r - (C6) 

J = d 2 g/dr u dr u l, u = r - (C7) 


Linear Measurement Model 

Given the large user/satellite separation and slowly changing geometry, com¬ 
mon practice is to employ an extended (linearized) Kalman filter that encompasses 
only the first two terms for g in Eq. (C3). The standard EKF equations for 
updating the a priori estimate of user position ( r u ) and clock bias (b u ) are as follows: 

*u = X~ + k(z ~ go) 
k = Ph T /(hP~h T + al) 

P + = (/ - kh)P {I - kh) T + k(vl)k T (C8) 


where 


and 


(C9) 

h = [(r s — r a ) T ig 1] (CIO) 



Also, P is the covariance of the a priori state estimate (x„), and <rl is the variance 
of the measurement error v, assumed to be Gaussian white noise. 


Nonlinear Measurement Model 

When PL measurements are introduced, the user/PL range is much less and 
more dynamic. Consequently, a Gaussian second-order filter 19 can be employed 
that accounts for the measurement nonlinearity with the quadratic component 
included 


z ~ go = hAx a + Ar£/Ar u /2 + v 


(Cll) 
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The corresponding GSO filter equations are expressed by the following: 20 
= x~ + k(z - go ~ r\) 

Ik = Ph T /(hP~k T + u 2 v + cr 2 ) 

P + = (I - kh)P~(I - khf + k(<j\, + al)k T 
-q = Tr[7P-]/2 

a 2 = Tr[JP~JP~]/2 (Cl 2) 

where x a , h , and J are as defined in this appendix. The new components account 
for the bias and measurement variance introduced by the quadratic nonlinearity. 
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